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Masses of heavy Standard Model fermions (top-quark, bottom-quark, and tau-lepton) 
play an important role in the analysis of theories beyond the SM. They serve as low- 
energy input and reduce the parameter space of such theories. In this paper Minimal 
supersymmetric extension of the SM is considered and two-loop relations between known 
SM values of fermion masses and running parameters of the MSSM are studied within the 
effective theory approach. Both 6-quark and r-lepton have the same quantum numbers 
with respect to SU (2) group and in the MSSM acquire their masses due to interactions 
with the same Higgs doublet. As a consequence, for large values of tan/3 parameter cor- 
responding Yukawa couplings also become large and together with tan [3 can significantly 
enhance radiative corrections. In the case of fe-quark two-loop O(ct^) contribution to the 
relation between running bottom-quark mass in QCD and MSSM is known in literature. 
This paper is devoted to calculation of the NNLO corrections proportional to Yukawa 
couplings. For the r-lepton obtained contribution can be considered as a good approxi- 
mation to the full two-loop result. For the 6-quark numerical analysis given in the paper 
shows that only the sum of strong and Yukawa corrections can play such a role. 
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1. Introduction 

One of the remarkable properties of the supersymmetric (SUSY) extensions of the 
Standard Model is the possibility to obtain nice unification of gauge couplings at 
the GUT scale. By means of one-loop renormalization group analysis of the Mini- 
mal Supersymmetric Standard Model (MSSM) in the beginning of 90s of the last 
century the scale of SUSY breaking (TeV) compatible with the unification at 10 16 
GeV was "predicted" ^. At present modern computer codes (SOFTSUSY^, SuSpect 
^, SPheno ^, ffmssmsc ^) routinely use two-loop renormalization group equations 
(RGE) to calculate the spectrum of superparticles given high energy input for SUSY 
breaking parameters. It is obvious that RGE at higher loops become a system of 
coupled differential equations so even to study gauge coupling unification one needs 
to know the value of other, e.g. Yukawa, couplings. 

Since for the moment the mass of an elementary particle is the only source of 
information about its coupling to Higgs boson (s), fermion masses are important 
low-energy input for all the models beyond the SM. 
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Top quark, bottom quark and tau-lepton are considered to be the heaviest 
fermions known in nature. In the context of the MSSM heavy (b and t) quark 
masses were studied in literature and leading two-loop relations between pole and 
running masses were found 6 7 8 . 

The pole mass does not depend on the renormalization scale /jl 9 l 10 lso the running 
masses can in principle be extractecQfrom it at any value of /2. However, in practice 
the scale should be tuned in a proper way to avoid large (high-order) radiative 
corrections. 

At the electroweak scale (~ Mz) the top-quark pole mass M t can be used to find 
the value of the running mass mP R defined in DR-renormalization schemdl 1 ! 12 ! 13 ! 
This is due to the fact that Mz ~ M t and there are no large logarithms in the 
relation. 

The story becomes more involved if one considers 6-quark and r-lepton. In a 
theory with many different mass scales m <C M it is not so easy to avoid the 
appearance of large contributions in the form of log M/p, and logm//2 with m,M 
corresponding to masses of light and heavy particles. This is a "non-decoupling 
feature" of minimal (MS-like) renormalization schemes ^3. In our case we have 
Mi,, M T <§; Mz and there are large logarithms in the relations. Moreover, for b- 
quark there is a renormalon ambiguity!^ that limits the precision of experimental 
pole mass determination. 

In order to solve the problem one usually employ the concept of effective field 
theory (see, e.g., Ref. [12] for a review) and perform "manual" decoupling of heavy 
particles. The procedure is well-known in the context of QCD Q2I ( see a i so a n j ce 
program RunDecPl). 

This approach allows one to relate running parameters in well-established ef- 
fective theory and corresponding parameters in a more fundamental theory by 
means of so-called decoupling constants which can be calculated order by order 
in perturbation theory. In some sense decoupling constants absorb leading contri- 
bution of heavy particles to various low-energy quantities. For the 6-quark two-loop 
0(a 2 s ) contribution to the relation between mjp R and mj^ s due to strong interac- 
tions was obtained within MSSM in Rcfs. [T!|J [SUl It is known from one-loop^ and 
two-loop^ calculations that strong corrections with virtual supersymmetric particles 
can be significantly reduced by contributions due to other interactions. In this paper 
the results for two-loop corrections proportional to Yukawa couplings of heavy SM 
fermions otf = y'j/^ir) with yf = {?/*,?/£>, ?/t} will be presented. 

Contrary to the quark masses leptons do not have large uncertainties due to 
confinement so tau-lepton pole mass can be extracted from the experiment with 
high precision M T = 1776.84 ± 0.17 MeV (see e.g. Ref. [22]). 

In principle this fact allows one to determine the value of the running mass (or 

a One should keep in mind that we are talking about experimental constraints on running param- 
eters of a model that allows one to reduce the parameter space, expressing, e.g., running masses 
of fermions in terms of of other parameters (heavy particle masses) 
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equivalently Yukawa coupling) very precisely. In the MSSM only one-loop super- 
symmetric contribution to the relation between pole M r and running m^? R masses 
is knowrPH. The value of the correction depends on parameters of the MSSM and 
in some cases can be of the order of 10%. 

Clearly, in comparison with the experimental uncertainty one-loop contribution 
is rather big. So it seems to be a good idea to calculate two-loop corrections. More- 
over, a general reasoning tells us that the inclusion of the two-loop result allows one 
to reduce the dependence of the final result on the renormalization (or decoupling) 
scale p,. 

Why decoupling constants are important in studying a theory beyond the SM? 
Together with renormalization group equations (RGE) they allow one to use a power 
of MS-like schemes in studying high-energy behavior of the MSSM. According to 
formal perturbation theory in order to obtain the value of, e.g., MSSM 6-quark 
running Yukawa coupling yi, at the GUT scale with L-loop precision one needs to 
perform a matching of an effective theory (e.g., SM or even Fermi theory) with 
more fundamental MSSM at (L — l)-loop level somewhere at the electroweak (or 
SUSY) scaleg So for one-loop RGE analysis decoupling constants are trivial and 
running parameters are continuous when one crosses a threshold of some heavy 
particle. The situation becomes more involved when two- or three-loop RGEs are 
employed. The parameters obtain a non-zero shift at the scale at which a heavy 
particle is decoupled. As it was mentioned above many codes use two-loop RGEs 
and one-loop decoupling corrections are incorporated. It should be pointed out that 
there exist a dilemma at which scale to decouple heavy particles and how many 
particles to decouple at chosen scale. 

The problem is that when we cross the threshold and decouple a particle we 
sometimes break a symmetry that guarantees the equality of coupling constants that 
enter different interaction vertices in a Lagrangian. For example, decoupling of only 
one squark breaks the supersymmetry which relates the interaction of quarks, gluino 
and squarks to that of quarks and gluons. As a consequence, one needs to introduce 
a new coupling constant in the effective theory without the squark. This coupling 
coincide with the strong coupling constant g s above the decoupling scale /idee but 
is not equal to g s below /idee- Of course, the difference can be calculated. However, 
when one goes from the MSSM to the SM and decouples every heavy superparticle 
at its mass a bunch of intermediate effective theories are produced with different 
symmetries broken (it can also be SU(2) symmetry of the SM) with different RGE 
equations and different threshold corrections. This way is certainly not the optimal 
one. In order to make use of all the symmetries presented in full theory during 
calculation of threshold corrections it seems to be a good idea to match the SM 
directly to the MSSM ( "common scale approach" ) ^ So we are left with the issue 
of choosing the decoupling scale. A common choice is Mz scale, (see, e.g. Refs. El 

b In principle, the result does not depend on the decoupling scale. Again due to truncation of the 
perturbative series one has to be careful when choosing a particular value. 
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I23[l . Clearly, at Mz high order terms can become important. Moreover, there exist 
some MSSM scenarios^ when masses of scalar superparticles significantly differ 
from that of fermion superparticles. In this case high order corrections can also 
improve the precision of calculation. 

In fact, in the context of MSSM RGEs are known up to three loopa^,. The anal- 
ysis presented in Ref. [25] was based on one-loop threshold (decoupling) corrections 
and it was mentioned the necessity of two-loop results for self-consistent study. The 
results presented in this paper together with that obtained earlieiGlEI] are aimed to 
partially fix this mismatch. 

The paper organized as follows. First of all, the approximation to the MSSM 
(so-called gauge-less limit) is described in Section [2j A special attention is paid 
to the tadpole diagram treatment in Section [3] Then a brief review (see Sec. [3} of 
decoupling procedure is presented . Section |5j is devoted to the results and numerical 
analysis of the calculated contribution in a wide range of parameter space of the 
MSSM. In the end of the paper Conclusions and Acknowledgments can be found. 

2. Gauge-less limit of the MSSM 

In order to simplify the decoupling procedure as an effective theory I considered a 
theory of free tau-lepton and five-flavor QCD with massive bottom-quark. Due to 
smallness of electroweak gauge coupling I neglected them. Moreover, the lightest 
Higgs boson is assumed to be much heavier then the bottom quark and tau-lepton 
so it is "left" in the MSSM. 

In the effective theory we employ MS minimal renormalization scheme and have 
the following set of running parameters: mass of the tau-leptorj^] (m^ s ), mass of 
the 6-quark (m^ s ) and strong coupling constant (ff^ 13 )- By means of well-known 
technique ^2 they can be related to the parameters of the MSSM defined in so-called 
DR-schemc0 

In such a simple setup decoupling corrections to the tau-lepton mass coincide 
with the correction^ that enter the relation between the and the pole mass 

M T . Due to strong interactions for the &-quark we have to take into account the 
difference between and Mf,. 

Since I neglected electroweak gauge interactions in the effective theory it is 
convenient to do the same thing in the MSSM. This approximation is called a 
gauge-less limit of the MSSIVPI 

In this limit there is no mixing between gauginos and higgsinos so only the 
latter have to be taken into account. Both charged and neutral higgsinos have the 
same mass that is equal to the absolute value of the supersymmetric Higgs mixing 

c In the considered approximation "running mass" coincide with the pole mass M T . 

d The issue of MS — > DR transition can be solved be decoupling of unphysical e-scalars in a way 

presented in Ref. 1191 or by a two-step procedure given in Ref. | 2 6l 

e Actually, with the first term of Large Mass ExpansiorEIHIlH] 
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parameter The mixing matrices for chargino (U and V) and neutralino (N) are 
given by the following expressions 
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It should be mentioned that in Higgs sector we can have a problem since in the 
MSSM quartic interaction of Higgs bosons are proportional to the electroweak gauge 
couplings. Nevertheless, in this paper I assume that there is a successful electroweak 
symmetry breaking so Higgs bosons have non-trivial vacuum expectation values and 
treat the gauge-less limit in a formal way. 

It is fair to say, that I am not going to be completely self-consistent within the 
limit. In the next section, Higgs sector of the gauge-less version of MSSM will be 
discussed together with an issue related to so-called tadpole diagrams. As it will be 
shown, in our case we have a very degenerate situation in Higgs sector with four 
massless bosons and four bosons with equal masses. Clearly, this is not satisfactory 
from the phenomenological point of view. 



3. Higgs sector and tadpoles 

In the MSSM like in the SM the electroweak symmetry is broken by a vacuum 
state, which can be characterized by the vacuum expectation values v\ and V2 of 
two Higgs doublets Hi and H 2 

In the true vacuum the first derivative of the effective potential (tadpole) should 
vanish. However, naive calculation of one-loop tadpoles for neutral CP-even Higgs 
bosons h and H shows that they do not vanish. This can be interpreted in the 
following way. The vacuum expectation values for the Higgs bosons that we substi- 
tuted in the tree-level Higgs potential minimize the tree-level potential instead of 
effective one. 

In order to deal with the problem I adopt the reasoning of Refs.|SJ|5T] Tree-level 
tadpoles for <fii and (f>2 that appear in the MSSM Lagrangian after substitution of 
j2) look like 

a 2 + a' 2 

Ti = (u 1 - v 2 )vi + mxVi - m 3 i> 2 , 

a 2 + n' 2 

rr, y T y , 2 2\ , 2 2 ro\ 

T 2 = (v 2 - v 1 )v 2 + m 2 v 2 - m 3 vi, (3) 



In this work I assumed that fj, > which corresponds the positive contribution to the muon 
anomalous magnetic moment. 
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where g, g' are gauge SU(2) x U(l) couplings (they will be neglected in what 
follows) and m\,m\, and m 2 are soft supersymmetry breaking parameters of the 
Higgs potentiaP. Both T x and T 2 are assumed to be non-zero and serve as counter- 
terms to cancel loop-induced tadpoles. As a consequence, I do not need to consider 
diagrams with tadpole insertions like in Refs. |3TJ |3H since they are precisely 
canceled by the tree-level counter-term. However, they do not disappear completely 
and manifest themselves in the mass matrices of all Higgs bosons. In the gauge-less 
limit mass matrices have the following form ($ = {<j)i,<j)f i = 1,2) 

( mi — TnT\ 9 /cot/3 — 1 \ , 



\ — 1 tan/3y 
SM *=[ V Ht 2 ) ( 5 ) 




Since T\ and T2 should cancel loop-induced tadpoles they are at least 0(g 2 ) 
with g being some coupling constant of the theory. Consequently, in order to obtain 
tree-level mass matrices one should set T± = Ti = and diagonalize only first term 
in Q. Clearly, this matrix has two eigenvalues (0, m§) and it is diagonalized by the 
rotation with angle f3. 

As it was mentioned in the end of the previous section this is not a satisfactory 
result. I decided to be slightly more close to the MSSM and introduce different 
masses for all physical Higgs bosons together with a different mixing angle a for 
CP-even states 



h J \(t>2/ V — sin 9 cos 



G+ 



= m ( if) • (8) 



After tree-level diagonalization linear and quadratic parts of the Higgs potential 
can be rewritten in the following way: 

V 2 = T H H + T h h + - (Mjj H 2 + Ml h 2 + Ml A 2 ) + M^ ± H+H~ 



bc+G- b G + H - \ ( G~ 
ba+H- b H + H -J \H~ 

Here H, h are neutral CP-even. Higgs bosons, A and G° — neutral CP-odd higgs 
and Goldstone boson correspondingly, H ± and G ± — charged higgs and goldstone 
bosons. Additional contributions to mass matrices come from SM$ (see Eq. (pi) 
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and can be expressed in terms of rotated tadpoles Th, Th 



T h J *< >UJ' (10) 
a vacuum expectation value v = \A>i + «|, and Higgs mixing angles (3, a ® 

(T H (cl,sp + s%cp) +T h s a c a s a ^p) , (11a) 

rS a -f3 + T h C a ^p) , (lib) 

,C a S a C a -[3 +T h {(? a Cp - S^Sp)) , (lie) 

- (T H c a -p - T h s a -p) , (lid) 

- (Ti? s a -p + T h c a -f3) , (lie) 
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(T H {slc a + cls a ) + T h (clc a - sls a )) , (llf) 

vs 2 p 

where eg = cos and sg = sin 9 for some angle 6. Since the contribution due tadpoles 
does not depend on external momenta it can be taken into account by introduction 
of additional non- minimal counter-terms for higgs masses and mixing. 

In the end of this section I would like to mention another issue related to Gold- 
stone bosons. In the gauge-less limit there are no gauge bosons to "eat" goldstones 
so the latter are massless and the global SU(2) x U(l) symmetry prevents them to 
acquire mass. However, in perturbation theory it is not so obvious. For example, 
if one does not take into account the contributions (111 during calculation of gold- 
stone boson self-energies T, GG (p 2 ) non-zero masses for G° and G are immediately 
generated since T, GG (0) ^ 0. It worth mentioning that we also run into a problem 
with spurious infra-red (IR) divergences when we try to calculate Feynman inte- 
grals with Y, GG (p) insertions. The role of tadpoles is crucial here since they precisely 
cancel this non-zero contribution keeping Goldstone bosons massless. 

In deriving the result presented here I was using non-zero masses for Goldstone 
bosons denoted by M Gq and M G± which allows me not to deal with mentioned 
IR problem explicitly. This can be justified since in the full MSSM when linear R^- 
gauge is employed Goldstone boson masses are proportional to that of corresponding 
gauge bosons. This introduces gauge dependence in the result. However, numerical 
analysis shows that the dependence is negligible if one takes into account tadpole 
contribution. 



4. Decoupling procedure 

The process of obtaining decoupling corrections for 6-quark mass is described in 
great details in Ref. |T5;. I closely follow the same procedure. Here I would like to 
stress some of important steps of the calculation. 
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In order to obtain a final result for the relation between running masses in 
the MSSM and our effective theory we need to know so-called "bare" decoupling- 
constants Cm f ,B up to two- loop level. These constants relate bare parameters in the 
effective five-flavor QCD (underlined parameters) to that of the MSSM 



nif.B = Cm/, b (a B , M B ) x m/,. 



f = {b,r} 



(12) 



and depend on bare coupling constants as = {a s , a*, at, Q t }b and on bare masses 
of heavy particles denoted collectively by Mb ■ 

The needed relations can be derived from ( 12 1 by proper renormalization of the 
left- and right-hand sides. Since TO/(/2) = rrij(p,) and rrifQi) = m^ R (p,) are defined 
in different theories we have 

m f,B - Z mj \a J x m f , 
m f , B = 

As a consequence an implicit equation 
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, DR. 



= {a s , at, at, a T } 



DR 



(13) 
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(14) 



can be solved in perturbation theory by expressing a in terms of a DR and ex- 
panding the result in a DR . 

Simple power counting tells us that in our case we need to know Z® R and Z^p~ up 
to one-loop level and both Z m/ 's - up to two loops. Renormalization constants for 
the MSSM parameters Z® R and Zffi can be found in Appendix A The expressions 
for Z^ that were used in my calculation looks like 
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(15) 
(16) 



with Cf = 4/3 and Ca — 3 being quadratic casimirs of SU(3) group. The renor- 
malization constants can be easily obtained, e.g., from the formula mj — yfV\/^/2 
(/ = 6, r) and well-known RGEs for Yukawa couplings yb, y T , and vacuum expec- 
tation value v\ (see Refs . 134) . l35 |) . 

Corresponding expressions for Z™? defined in five-flavor QCD can be found in 
Refs. O [Ml As it was mentioned earlier Z^ s = 1 in our approximation. 
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Since there is a^ s in the right-hand side of ( 14 1 we also need one-loop relation- 
between a^ s and a^ R . This relation together with the expressions for full MSSM 
one-loop threshold corrections can be found in Appendix D of Ref. [2TJ In order to 
derive needed contribution to the fermion mass decoupling constants self- 

energies presented in this reference have to be expanded in external momenta and 
masses of the considered fermion and the mixing Q has to be taken into account. 

We are left with the only missing piece of the formula (14 1, i.e., two-loop contri- 
bution S(m] to Cm/ - Corrections of the order 0{a1) were found in Refs. [19| [20l The 
contribution proportional to Yukawa couplings of heavy SM fcrmions is obtained 
in this paper by means of a FORM 37 program specially written for calculation of 
decoupling constants in the MSSM. Corresponding diagrams were generated with 
the help of FeynArtsCHH 



0{cfi,5<£\Mz),% 0(a s a f + a}),5(l 2) (M z ),% 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

mo, TeV mo, TeV 



Fig. 1. Two-loop 0(a 2 s ) and 0(a s aj + a 2 ) contributions to &(,m\ as functions of mo and mi/% 
for Aq = 0. Upper row is for tan j3 = 10 and the lower one is for tan j3 = 50. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

mo, TeV too, TeV 



Fig. 2. Comparison of one- and two- loop contribution to the decoupling constant for the fe-quark 
running mass ( mb as functions of rag and nij^ for Ag = 0. Upper row is for tan/3 = 10 and the 
lower one is for tan /3 = 50. 

5. Results 

In this section I present the numerical analysis of the obtained result. Analytical 
expressions are huge and have been stored in the form of Mathematica code and 
GiNaC archive formal 8 The former allows to obtain numerical value for 6 (ml and 

(2) — 

given the running parameters of the MSSM. The latter gives us an opportunity 
to include the calculated corrections in the SoftSusy program ^ to evaluate them 
and to see the influence of the two-loop thresholds on the resulting spectrum and 
running parameters. In SoftSusy matching performed at the electroweak scale so 
for numerical results it is assumed that decoupling scale /id cc = Mz- 

In Fig. fTlone can find typical dependence of different two-loop contributions to 

— (2) 

the decoupling constant S( m l of the fr-quark mass on Constrained MSSM parameters 



s Both the expressions are available from the author by demand. 
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Fig. 3. Comparison of one- and two-loop contribution to r-lepton mass decoupling constant 
as functions of mg and m^/ 2 for Aq = and tan/3 = 50. 



mo and m,\ /2 for fixed value of ^4o = and for tan /3 = 10 (upper row) and tan /3 = 50 
(lower row). 

As one can see corrections due to Yukawa interactions tend to compensate 0{a s ) 
contribution. The effect of 0(a s otf + ocj) increases with tan/3. The comparison of 
first column of Fig. [I] and the resulting two-loop corrections presented in Fig. [2]) 
shows us that for large tan (3 — 50 total S£m b varies in the range of 2 - 4 % while 
0(a 2 s ) contribution varies in considerable wider range 2-11 %. From this fact one 
can immediately deduce the importance of two-loop Yukawa decoupling corrections 
for 6-quark running mass in the region of large tan (3. 

The sum of the above contributions is compared with full one-loop MSSM decou- 
pling constants in Fig. [5J It is easy to see that the resulting two-loop correction both 
for large and low values of tan (3 lies in the region of few percents and does not ex- 
ceed the relative uncertainty of the input parameter fhb = m^ IS (m^ ls ) = 4.20±0.18 
GeVEa 

In some sense it is a bad news since we obtained the result that is negligible. 
Nevertheless, one can be sure that large one-loop threshold corrections for 6-quark 
running mass widely discussed in literature (see, e.g., Ref. I4ip are indeed reliably 
approximate the full result. 

For tau-lepton strong interactions do not contribute to the mass decoupling 
constant at two-loop level so the corrections are relatively small in this case. Typical 
value of two-loop corrections 0{a 2 + a T a q ) with q = t,b for low values of tan/3 is 
10~ 3 % which is small even with respect to relative experimental error of the pole 
mass M T (10 -2 %). With the increase of tan (3 considered contribution is enhanced 
and can reach the value of few percents (see Fig. [3| 

In the end of this section I would like to demonstrate the role of tadpole contribu- 
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_4 I 1 1 1 1 .15 I 1 1 1 1 

500 1000 1500 2000 500 1000 1500 2000 

M Go , GcV M Go , GcV 

Fig. 4. The dependence of two-loop 0(a s af + ay) correction to ^Cm^ on Goldstone boson mass 
M^ Q for SPSla (mo = 100 GeV , m 1/2 = 250 GeV , A = -100 GcV , tan/3 = 10) and SPS4 
(mo = 400 GeV , mirt = 300 GeV , Aq = 0, and tan/3 = 50) scenarios. Contribution due to 
tadpoles (dot-dashed line) almost completely cancel (dashed-line) naively calculated correction 
without tadpoles. 

tion discussed in Sec. [3] Figure jijshows the dependence of 0{a s af+a 2 ) contribution 

to 8(ffl b on Goldstone boson mass Mq for SPS1 and SPS4 scenarios^. One can see 
that if we neglect tadpole contribution we overestimate the value of the correction 
and introduce significant dependence on Mq which can be interpreted as a gauge 
dependence. 



6. Conclusions 

The biggest experimental facility in the world, Large Hadron Collider™, has already 
been built and we are waiting for first physical run of the machine. We hope that 
there will be something that allow us to solve at least some of the problems of the 
SM and we believe that it will be supersymmetry. 

The MSSM is a viable candidate for a theory beyond the SM. It has a lot of 
parameters most of which related to supersymmetry breaking and they have to be 
determined from future experiments. However, other parameters are already con- 
strained from known low-energy input. Both bottom-quark and tau-lepton masses 
are among them. 

In this paper the relations between corresponding running masses m]p R and 
to^? r defined in the MSSM and known low-energy experimental input (mt and M T ) 
were considered. Two-loop (decoupling) corrections to this relations proportional to 
Yukawa couplings of heavy SM fermions were calculated. 

For the b-quark it was found that 0(a s ctf + a*) corrections are important for 
large tan f3 and significantly reduce two-loop strong contribution calculated earlier 
l 19 l 20 [ From Figure [l] one can see that mentioned contributions usually have different 
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signs. With the increase of m\ji (mo) absolute value of the corrections increases 
(decreases). Due to this kind of behavior they tend to compensate each other in 
the whole region of considered too — plane lowering the sum below current 

uncertainty in the input parameter fhb- 

Two-loop threshold corrections to r-lepton mass that were obtained in this paper 
turns out to be negligible in the region of low tan (3. With the increase of tan /3 they 
can reach the value of few percents and exceed the experimental error of M T . 

From the presented analysis it is obvious that the two-loop decoupling correc- 
tions are too small and can not significantly modify MSSM spectrum produced by 
public computer codes. Variations in the spectrum due to calculated corrections are 
comparable with variations due to uncertainties in the low-energy input parameters. 

Nevertheless, obtained result allows us to be sure that, e.g. for large tan/3 one- 
loop approximation is good enough. Moreover, if some new physics is established at 
LHC we will be ready to perform precision tests of SUSY models and their GUT 
extensions with the help of three- loop RGEs. 
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Appendix A. Renormalization constants 

To obtain finite result for threshold corrections we need to rewrite bare decoupling 



constants (14) in terms of renormalized parameters of the MSSM. I collected all 
the needed counter-terms in this Appendix. In what follows the following notations 
are used: mg corresponds to gluino mass, Af with / = {t, 6,r} — soft trilincar 
couplings, X t — A t — /i cot /3, and Xf = Af — /i tan (3 with / = {b, t} are off-diagonal 
elements of sfcrmion mixing matrices. Sfermion masses arc denoted by m% and 

/l,2 

to? (sneutrino). For sines and cosines of sfermion 9 f mixing angles abbreviations 
Sf = sin Of and s n f = sin(n6*y) (n = 2,4), etc are used. 

For bare scalar mass m 2 B , bare fermion mass Mg, bare coupling constant ag, 
and some bare mixing angle 9g corresponding counter-terms 8m 2 , SZm, 6Z a , and 
69 are defined by m 2 B — m 2 + 5m 2 , Mb = (1 + SZm)M, as = a(l + SZ a ), and 
9b = 9 + 66 with to 2 , M, a, and 9 being running parameters in DR scheme. 
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Appendix A.l. Scalar quarks 



(47r) 5m 2 ^ = C F y [i m i 2 _ TO g 1 ) s 2fc - 4 ( m l + m l- s 2b m b m g )j + -ym b X T s 2b 
+ yc 2 { A t + m l + m l s t + m l c t + M H ±c l + M G±4) 

+ { A b + m l + m i c t + m i s * + M "±4 + K±4) 

+ °J ( m b c b - m t s t + 3s 2 bin b A b + m^(l + 2s 2 2b ) + m? 2 (l - 2sl b )^ 

+ ^ (mlsl m 2 c 2 + m b .s 2b A 4 ) - ^ (% 2 + % 2 ) , (A.l) 



(47r) <5to? 2 = C F y (( m ^ 2 _ m f 1 ) s 26 _ 4 ( m 6 + m | + s 2b m b m g )j - -^m b X T s 2b 

+ ( A * + m i + m ^ s * + m l£ + M M + M G±4) 

+ yc 2 (a 2 . + m? 2 + m\c 2 t + m| 2 s t 2 + M 2 H± s% + M 2 ±4) 

+ y [a\- M 2 + ^ (M 2 S 2 + M^cl + M 2 A s 2 + M 2 „ C 2 )) 

+ ^ [mlsl - m 2 t c 2 b - 3s 2b m b A b + m|(l - 2s 2 2b ) + m \{\ + 2«§ 6 )) 

+ y ( m fc c 6 _ m t s 6 _ m b s 2b A t ) - ,u 2 ^Y s b + °J c b) . ( A - 2 ) 



( 4 tt) <W b 



- (A 



C' F — ({rnl 2 - m ~ bl )-Y + ^m b m g c 2b ^j + —m b X T c 2b 

(A 2 - M 2 + m 2 / t + m\/ t + M| ± s 2 + M 2 ±( 

((ro^ - m? 2 )s 4h + 3m b A b c 2b ^ + ^-m b c 2b A 
(a b a t \ ( 

VT'TJK 



e 2 
Ctb s 2b 

TIT 

Ob_ 

e 

s 2 b 
4 



+ J7i| 2 



2™; 



9 9 
»l b 2 



Til 



m 



(A.3) 



The counter-terms satisfy the following relation 



SMq = 8m^c 2 + 8m 2 ^s\ + (mf 2 — m^)s 2b 59 b — 2m b 5m b 

= Sm-^c 2 + Smf^s 2 + (m 2 ^ — m 2 ^)s 2t S9t — 2m t Sm t , (A.4) 
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which is a consequence of SU (2)-invariance of the MSSM Lagrangian. By means of 

2m b X b 



S2b 



C2b 



bi 



Ml 

b 



(A.5) 



it is easy to convince oneself that (A.l) (A. 2) and (A. 3 1 can be rewritten in the 
following way 



(4n)50 b 



wherc0 



Srrib 
m b 

sx b 

X b 



Mi 



SMj 



SMi ± 





Ml 

b 



2 

mf 

bi 



9 

mf 

b 2 



SMi 



Am 



, Srrih 



± 



2m b X b 



b "^2 

m b mf 



(Sm b X b + m b SX b ) 



bi 



m b X b 



/ 5m b 

V m b 



bi 

sx b 
x b 



(A.6) 
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lot 
2~e 



2C F 
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Xh ( 2 
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e 
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e 
3 at 
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.4 



<5M? =-4C F — m? + 



«6 

e 
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M- 
Q 

■2 
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M% 



a i 



Mi 



M~f 



A T - 



b 



£ 



£ 



A/ 2 - 

Q 



M? 



SMi 



AC F —m, 
e 



A 



A; 



-Hi 



A: 



SMi 

b 



(A.7) 

(A.8) 
(A.9) 
(A.10) 
(A.ll) 
(A.12) 
(A.13) 

(A. 14) 
(A.15) 



Renormalization constants for top-squark masses and mixing can be obtained from 
the expressions above by substitution b <-» t, cp <-> sp, c a <™> s a and a T — ► 



h needed counter-terms can be extracted from corresponding beta-functions given, e.g., in Ref. 1451 
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Appendix A. 2. Scalar leptons 



(47r)<5m? ^ a ^(A 2 T -n 2 + m 2 fi s 2 T + m\c 2 T + m% - m\ + c 2 M 2 ± + s 2 M 2 H± ) 



(A.16) 



(47r) 5m 



- -^s 2 T (A 2 T - p? + m\ + m? + M 2 ± c 2 + M 2 H± o fj 

£ c 



' 2 sV 



3X b m T s 2r 



+ — (m 2 fl (l + s 2 T ) + m? 2 (l - a\ T ) + 3s 2r A T m T + m 2 r c 2 T ) (A.17) 



(47r) dm 2 



a 



T c 2 T (A 2 . -p 2 + m\ + m% + M 2 ± s} + M 2 H± c 2 Si ) - —3X b m T s 2T 



" 7 



+ ^ (m 2 fi (l - s 2 T ) + m? 2 (l + a| T ) - 3s 2r A T m r + m 2 T s 2 T ) 



(A.18) 



(4tt) <5<9 T = 



3— m T X b c 2 i 

e 



+ ^ S -f(A 
e 2 



S2t a T 
4 e 

2 i ™2 



(m^ + m\ - 2m 2 ) 
/r + ml + M 2 H± sl-< nf2 - 
+ ^ ((m^ - m? 2 )s 4T + 3m r A r c 2 r) 



1 

m\ - m\ ' 



(A.19) 



Appendix A. 3. Higgs sector 

In order to obtain renormalization constants for Higgs masses and mixing angles 
we need to consider divergent parts of self-energy diagrams. The straightforward 
calculation leads to the following results 

(4tt) 5m% i = 3%g (m 2 h + m\ - 2m 2 + m| 2 + (A t - ^ ) 2 
+ 3^c 2 (ml + ml - 2m 2 + m| 2 + (A b - ^ ) 2 
+ ^c 2 (m\ + m\ 2m 2 T + m| 2 + (A b ^ f) ( A. 20) 

£ \ CQ 



(4tt) Sm% 2 = 3 — c 2 (ml + m 2 ^ - 2m 2 + + (A t + fi— )'■ 
£ y cq 

+ (ml + m? 2 - 2m^ + mj, + (A + ^ ) 2 ^) 

+ ^s 2 (m?^ + m? 2 - 2m 2 + r< + (A 6 + ^ ) 2 ) (A.21) 
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(4tt) 59 



,Oib 



(ml + m? - 2m\ + (A t - M - ){A t + ^ ) 
\ 1 2 s e c e 



se 



ce 



3— Im + mf -2mi + (A b -[i— )(A b + ^ — ) 



Q r 



2m2 + (A T - M ^)(A r + M ^) 
eg se 



, a i 



_ a. 



cese 



(A.22) 



Here <&i = {G,G ± ,H}, <f> 2 = {A H ± ,h}, 9 = {/3, (3, a} correspondingly. It is easy 
to see that for 5(5 one obtains different expressions when considers CP-odd (<5/3o) and 
charged higgses (5j3±). This apparent problem "is solved" if we remember that in 
the true gauge-less limit M^ q — M^ ± — and M\ = M'j I± so 5f3 — 5(3±. However, 
it is convenient to keep the masses different and instead of renormalization of bare 
angles /3q and j3± use the following explicit counter-terms 



8CZ 



{M% - M 2 A ) 50o x G A Q 
(M% ± - M% ± ) 5P± x (H+G- 



■G+H-), 



(A.23) 



where 5[3 and 5(i± are given by (A.22 1 for (<E>i, *E> 2 ) = (G,A ) and (4>i,<& 2 ) = 
(G ± ,H ± ) correspondingly. I would like to mention that it is not the end of the 
story. For the moment we neglected the contributions due to tadpoles described in 
Sec! 



(4tt) Stm^ = 



ml + ml - 2m? + (A t - ^){A t - 4") 
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C 
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C 



Sg Sp 



(A.24) 



(A.25) 
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(4tt) 8 t 



3^ (ml + ml - 2m\ + (A t ^ ){A t + [i— ) 
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It is interesting to note that for $i = {G, G*} and $ 2 = {A, i? ± } 
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(A.29) 



where for m% 



l G,G± 



(true goldstoncs) and m% 



m. 



m 



H± 



m\/(cpSp) = B^/{cpSfj) the result can be cross-checked with the help of one-loop 
renormalization constants 



3 aft 1 a T 
~2~ ~ 2~£~' 



(4ir)SB = 3A, 



3 A, 



(47r)<5f2 



3 at 
'2 £ ' 



and relations 



Sf3 = cos 2 /3 S tan j3 



(A.30) 
(A.31) 

(A.32) 



Appendix A. 4. Other renormalization constants 

For completeness I present one-loop counter-terms for couplings 



(4n)6Z a . = (6 - 3C A ) 



(4ir)5Z ab 
(4Tr)6Z aT 



-4Cj 

-4C F 
, a t , 
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(A.33) 
(A.34) 
(A.35) 
(A.36) 
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and for the masses of i-quark mt and higgsinos m x (6-quark and r-lepton are con- 
sidered in Sec. [3} 



(4tt)(5Z„ 

(4TT)SZ rn 



-2C F - 



3 at 1 at 

2Y + 2T'- 



3 fat, 
2 V e 



(?+?) 
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(A.37) 
(A.38) 



Clearly, 5Z m coincide with the renormalization constant for /i (A.ll). 
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